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1  Introduction 


The  investigation  of  polymers  is  one  of  the  fast  developing  field  of  today’s  mate¬ 
rial  science.  Besides  the  large  variety  of  experimental  methods,  the  theoretical 
desription  of  the  polymers  becomes  more  and  more  important.  Earlier  only 
semiempirical  methods  of  quantum  chemistry  were  applicable  to  the  polymers 
with  larger  unit  cell,  but  recently  with  the  very  fast  development  of  computers 
the  first  principle  (ab  initio)  methods  are  also  used  for  these  systems.  The 
most  widely  applied  ab  initio  method  to  investigate  periodic  polymers  is  the 
Hartree-Fock  crystal  orbital  (HF  CO)  method  has  been  developed  initially  by 
Lowdin  [1],  and  further  introducing  an  LCAO  expansion  by  Del  Re,  Ladik 
and  Biczo  [2]  and  Andre,  Gouverneur  and  Leroy  [3].  Several  HF  CO  program 
packages  are  available  today.  We  mention  here  the  Crystal  series  of  programs 
from  Pisani,  Dovesi  and  Roetti  [4]  and  the  program  developed  in  Erlangen  by 
P.  Otto  [5]. 

Often  the  fundamental  gap  of  a  polymer  is  of  interest  since  this  determines 
the  conduction  properties  of  the  system.  It  has  long  been  realized  that  at  the 
HF  level  the  band  gaps  of  polymers  are  overestimated  and  thus  considerable 
effort  has  been  concentrated  on  developing  techniques  to  ccdculate  electron 
correlation  effects.  In  infinite  systems  like  polymers  a  method  must  be  chosen 
which  is  size  consistent  and  therefore  many  body  perturbation  (MP2,  MP3) 
[6,  7,  8]  and  coupled  cluster  techniques  [9,  10]  have  been  applied  to  treat  these 
systems.  Calculation  of  the  quasi  particle  (QP)  band  structures  of  polymers 
using  the  inverse  Dyson  equation  [11]  allows  a  more  accurate  prediction  of 
band  gaps. 

A  new  implementation  of  QP  method  described  in  the  first  Report  and 
pulished  in  J.  Comput.  Chem.  makes  these  calculations  much  faster  [12]. 

In  the  second  Report  we  presented  the  first  version  of  the  direct  SCF  cristal 
orbital  program,  described  some  possible  solutions  of  the  convergency  prob¬ 
lems  which  arise  in  the  SCF  calculalations,  and  reported  the  first  results  of 
correlated  band  structure  calculations  for  polymers  with  leirge  unit  cells. 

The  most  time  consuming  part  of  the  quasi  particle  band  structure  calcula¬ 
tions  is  the  integral  transformation  (the  so-called  four  index  transformation). 
In  this  report  we  present  an  algorithm  that  reduces  the  required  time.  Our 
program  packace  makes  it  possible  to  calculate  the  correlated  band  structure 
of  polymers  with  large  unit  cells.  We  tried  to  find  different  fields  of  polymer 
science  from  synthetic  metals  to  biological  systems  to  demonstrate  the  wide 
range  of  application  of  our  program. 

The  discovery  of  electrically  highly-conducting  polymers  is  one  of  the  signif¬ 
icant  developments  of  recent  times.  These  materials,  also  called  synthetic  met¬ 
als,  combine  the  electric  properties  of  metcils  with  the  advantages  of  polymers 
such  as  smaller  weight,  greater  workability,  resistance  to  corrosion  and  lower 
cost.  The  hydrocarbons  with  ladder  structure  play  very  important  role  among 
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these  materials,  one  representant  of  these  family  is  the  poly-(pert-naphtalene) 
(PPN)  .  The  manufacture  of  light  emitting  diodes  based  on  conducting  poly¬ 
mers  has  been  a  breakthrough  in  polymer  electronics.  For  optical  applications 
the  most  promising  materials  are  among  others  poly(p-phenylene  vinylene) 
(PPV)  ajid  the  poly-(p-phenylene)  (PPP).  DNA  plays  an  important  role  in 
living  systems  carrying  the  genetic  information.  The  electronic  structure  of 
DNA  is  also  very  important  from  the  point  of  carcinogenesis  [13].  In  the  last 
part  of  the  Report  we  describe  the  results  of  calculations  for  nucleotide  base 
stacks  (in  B-DNA  structure),  PPP  ((Cgffe)*)  and  PPN 


2  Theoretical  Background 

2.1  Quasi  Particle  Band  Stucture 

The  procedure  for  the  calculation  of  the  quictsi  particle  band  structures  of 
polymers  has  been  developed  over  several  years  in  Erlangen,  initially  by  Suhai 
[6]  and  then  by  Liegener  [7]  (see  eJso:  [8])  as  a  generalization  of  the  electron 
polaron  model  [14].  The  introduction  of  Green’s  function  formalism  allowed 
a  more  accurate  determination  of  the  quasi  particle  energies  by  virtue  of  the 
Dyson  equation; 

G(cj)  =  Go(u;)  +  Go(tc;)S(cc^)G(a;) 

G(ct;)  is  the  Green’s  matrix  for  the  perturbed  system  and  Go(ci;)  is  the  corre¬ 
sponding  Green’s  matrix  for  the  unperturbed  system,  which  is  given  by 

Go(a;)  =  (a;  •  1  —  £)  \  (1) 

Here  5  is  a  diagonal  matrix  containing  the  HF  band  energies,  uj  is  the  quasi 
particle  energy,  S  is  the  self  energy  term.  Using  only  the  diagonal  elements 
of  the  second  order  self  energy  matrix,  the  quasi-particle  band  energy  of  state 
I  satisfies  the  scalar  inverse  Dyson  equation  [15] 

Luj  ^  ei  +  (2) 

where  uji  is  the  quasi-particle  band  energy,  ei  the  corresponding  HF  band 
energy  and  the  diagonal  second-order  self  energy  term.  It  should  be  noted 
that  in  the  polymer  case  the  above  equations  are  written  in  the  CO  basis  and 
thus  the  index  /  is  a  composite  index  and  refers  to  both  a  band  index  and 
k-index^  ie,  I  =  (iyki).  The  expression  for  the  self  energy  term  comes  from 
many  body  perturbation  theory  and  can  be  written  expanding  the  composite 
indices  in  the  CO  basis  as, 


X 


iCocc  k^ 

k<l^occ 

2  *  d"  2  •  2  *  Re[u,*j7.^]  •  Re[‘U|j7/j]  2  *  Im.['i7i’jfc/]  *  Irn[‘Uij7^] 


w(t,  ki)  +  £(;,  kj)  -  e{k,  kk)  -  e{l,  k)  +  irj 


h/  ^ici+kj,kif+ici 
jiocc  jtfc  ki 
k<l^occ 

2  •  +  2  •  \vijik^  —  2  •  •  Re[t;,-j7fc]  —  2  • 

u{i,  ki)  +  £(;,  kj)  -  e{k,  h)  -  e{l,  k)  -  i-q 


(3) 
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where  the  integral  Viju  is  defined  as 


The  factor  /«  =  1  if  fc  7^  1  and  fki  =  1/2  if  fc  =  i.  It  is  introduced  because 
of  the  restriction  of  the  summation  over  the  band  indices  k,l.  (f>i{k^fi)  is  a 
solution  of  the  polymeric  HP  equations  expressed  as  a  linear  combination  of 
Bloch  orbitals, 

<i>i{k,f)  =  Y,Cia'>Pa{k,f)  (4) 

a 

and 

Mk,  n  =  E  -Ri-  R.),  (5) 

I 

where  {Xa(^~  Rl  —  ■Sa)}a=i,2 . N  is  a  set  of  basis  functions,  Ri  the  lattice 

vector,  Ra  is  the  center  of  a-th  basis  function. 

From  the  expressions  described  above  there  are  two  methods  which  may 
be  adopted  for  locating  the  quasi  particle  energies.  The  first  method  searches 
a  grid  of  energies  for  graphical  solutions  of  the  inverse  Dyson  equation  as  was 
proposed  by  Liegener  [7],  The  method  of  Palmer  and  Ladik  [12]  locates  the 
solutions  by  iterating  until  a  convergence  criterion  is  met.  Using  equation  (3) 
the  self  energy  may  be  calculated  for  the  band  and  k-point  of  interest  by  setting 
the  quasi  particle  energy  equal  to  the  HP  energy  in  the  zeroth  iteration.  The 
result  is  then  used  in  the  inverse  Dyson  equation  (2)  to  calculate  the  new  quasi 
particle  energy  and  the  process  is  repeated  until  convergence  is  reached.  For 
the  computational  details  see  [12]. 

2.2  Integral  transformation 

In  a  large  part  of  the  molecular  quantum  chemical  calculations  of  correlation 
energy  the  bottleneck  is  the  integral  transformation  from  the  original  (mainly 
Gaussian)  basis  set  to  the  MO  basis.  For  polymers  this  problem  is  more 
serious  bacause  of  the  ^-dependence  of  the  CO  coefficients.  We  describe  here 
the  polymer  analogue  of  the  molecular  transformation  scheme. 

2.2.1  Direct  transformation 

The  calculation  of  the  self  energy  term  requires  the  transformation  of  the 
two-electron  integrals  in  the  atomic-orbital  (AO)  basis  to  the  CO  basis  -  in 
molecular  case  this  is  the  so-called  four  index  transformation. 

ix°y^\x"xi)  (6) 
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The  CO-s  axe  expressed  in  the  AO  basis  in  the  following  form: 

a  I 

We  can  get  the  rough  form  of  the  direct  transformation  formulas  by  substitut¬ 
ing  this  expression  into  the  right  hand  side  of  (3). 


a/3^6  UHL 


To  transform  all  the  CO  two-electron  integrals  the  number  of  operations  is 
proportional  to  (n/cp^  *  726/^  *  ncell^  *  nbf^).  Here  nbf  is  the  number  of  basis 
function  used  in  (5),  ncell  is  the  number  of  interacting  cells,  nkp  the  number 
of  k  points  used  in  the  numerical  integration. 


2.2.2  Transformation  through  Wannier  functions 

Suhai  [6]  proposed  an  intermediate  transformation  step  to  solve  this  problem. 
The  AO  integrals  are  first  transformed  to  a  Wannier  function  (WF)[16]  basis 
and  in  the  second  step  the  WF-s  are  transformed  to  the  CO  basis. 


(xWlx^Xi)  — ^ 


ty^-s  are  the  WF-s  constructed  from  the  COs  by  the  Fourier  transform: 


k 


(7) 


Using  the  fact  that  besides  all  functions  are  also  eigenfunctions  of 

the  HF  Hamiltonian,  we  have  a  possibility  on  one  hand  to  create  real  Wannier 
functions  (with  the  choice  A“^  —  ~^n)  other  hand  to  reduce  the 

spatial  extension  of  the  WF-s  [6,  8].  In  our  calculation  we  used  the  method 
of  Knab,  Forner  and  Ladik  [17]  that  chooses  A*-s  in  a  way  to  maximize  the 
WF  norm  contributions  from  the  0-th  cell.  Using  (5)  and  (4)  the  Wannier 
functions  can  be  expanded  in  the  set  of  xLi^ ~  ~  -^a)  basis  functions: 


“■t  =  E  E  -Rj-R.),  (8) 

a  J 


where 


(u,?u,/Kwf)  =  E  E  4-“4-'4-"4-‘(xixiix?xl).  (9) 

a/S-ri  UHL 
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Here  /,  J  etc.  are  running  indices  and  J,  H  etc.  are  fixed  cell  indices  into 
which  the  different  Wannier  functions  are  localized.  The  number  of  required 
operations  is  now  proportional  to  [ncell^  *  nbf^  *  ncell^  *  nbf*) 

In  the  second  step  we  transform  the  WF-s  to  CO-s  back: 

=  N-i 

L 

Using  this  transformation  we  can  calculate  the  CO  integrals: 

A.B,C 

The  number  of  required  operations  is  proportional  to  {nkp^  ^ncell^  ^nbf^). 
We  can  see  comparing  the  number  of  required  operations  this  and  in  the  direct 
transformation  case,  that  it  is  worth  to  use  the  above  given  method  if  the 
number  of  cells  taken  into  account  is  smaller  then  the  number  of  k  points 
which  is  the  case  in  our  calculations. 

2.2.3  AO  to  WF  integral  transformation 

In  the  HF  CO  program  [5]  the  atomic  twoelectron  integrals  are  calculated 
only  for  a  nonredundant  set  of  cell  quadruplets  (  OJHL  ).  It  means  that  the 
cell  indices  are  restricted  to  a  certain  neighbour  of  the  reference  cell  and  the 
trivial  symmetry  of  the  two-electron  integrals  is  taken  into  account  (for  the 
details  see  [18]).  We  describe  here  the  form  of  integral  transformation  formula 
(9)  using  this  non-redundant  set. 


“/S-r-J  IJhL 


Hs 


n  jo  ka  10 
jH-QjL-JjJ-Hj-L 


^k0 


^la 


tL-0  lii-J  iI-H  iJ-L 


^ka  ^10 


+ 

+ 

+ 

+ 

+ 

+ 

+ 

X  i  •  (xixilxfxi), 


where  the  *  shows  the  restricted  summation  for  the  non-redundant  integrals 
and  t  =  (1  —  \Sjj5ai3){'i-  —  ~  In  the  next  step  we 
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use  the  translational  invariance  of  the  (x{(l)x^(l)|xf^(2)x|'(2))  two  electron 
integrals; 

(xi(l)xi(l)lxf(2)xh2))  =  (x”(l)xi-'(l)lxf'(2)xi-'(2)). 

Using  the  J'  =  J  —  I,  H'  =  H  —  I  ...  notation  we  obtain: 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

X  ^  ■  (XaX^lxf  xf )  =. 

where  the  *  means  the  same  as  before.  After  some  rearangement: 


E  E 

*{E*( 

X 

a/37^  J'H^L’ 

I 

( 

+ 

( 

X 

( 

xf-  (XaX^lxrxf) 

a/37<y  J'H'L^  ! 


ja 


‘■ky 


i-OjJ'-J+I  jL’-H+I  jH'-L+I 


^kS 


^^.o,Id^-J^!dl-Hdr-L^! 

^^-o^!d^-j^!di--«^!dt^)} 


From  a  computational  point  of  view  it  is  worth  to  do  this  transformation 
in  four  separate  steps.  In  each  step  we  transform  only  one  AO  to  a  WF. 

•  1st  step: 


(x°xi'Ixft^? 

(x°xflxf4 

(x°xi'lxf^.^ 

(x°xflxf^f 


)  =  E  E*E*^  •  4  “■^^(xaxi'lxfxf ), 

5  L'  ! 

)=ErE‘‘-'i?;-^"'(xWixfxf'), 

6  L'  J 

)  =  EE*E*i-4-^"'(x°xi'ixfxf), 

S  L'  i 

)  =  ErE’i-4-""'(x°xflxfxf). 


6  U  J 


r 
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•  2nd  step: 


(x°xi>N/)  =  EE’E*{  <-^"'(x°x,^'lx?V)  + 

7  / 


(x°xi>f-,")  =  EE  E  {  <-""^(x°x,^'lxf-f)  + 

7  H'  j 

dky~^^^ixlxi'\x?'wf')} 


•  3rd  step: 


!-?»/)  =  EET'<ii:»-''*'(x2xiV?»/) 

P  r  1 


(xWlt.>/)  =  EE'E‘<'"^'(x2x^>?»/) 

fi  J‘  1 


(x»?^h  =  EE‘E'<'r°"'(x2xi'K»f) 

/3  r  ! 


(x^/l^f-f)  =  EET‘4-""^(x“x^'K<»f) 

13  J'  i 


•  4th  step: 


KM  hW)  =  E  E*{ 

“  I 


4a^  (Xa^fk.M)  + 

4;^  (x>fkW)  + 
4'"  (x>°kW)  + 
<^ia  (xX^kW)} 


The  number  of  required  operations  is  proportional  to  ~  4*{nbf^  *ncell^), 
which  is  ncell^  times  more  than  in  monomer  (molecular)  case. 
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3  Results 


With  the  extremely  fast  development  in  computer  science  it  is  possible  to-day 
to  taylor  polymers  with  different  required  properties.  Our  goal  is  to  develop 
a  program  package  that  is  a  proper  tool  for  this  purpose.  In  this  section  we 
report  the  results  of  the  program  developments  during  the  last  year  and  the 
computational  results  for  different  polymeric  systems. 

3.1  Changes  in  the  integral  transformation  program 

In  the  latest  version  of  the  QP  program  the  computational  time  problem  of  the 
solution  of  the  inverse  Dyson  equation  was  solved  using  an  itterativ  method 
[12].  The  applicability  of  this  QP  program  was  limited  by  the  I/O  and  time 
requirements  of  the  AO  to  WF  integral  transformation.  This  section  we  sum¬ 
marize  the  changes  that  solve  these  problems. 

•  .4s  we  have  seen  in  section  2.2.2  from  computational  point  of  view  it 
is  worth  to  use  Wannier  functions  in  the  integral  transformation.  The 
quality  of  WF-s  is  very  important  in  this  process.  The  smaller  is  the 
coefficients  in  the  subsequent  units  as  compared  to  the  reference  one, 
the  smaller  is  the  loss  in  the  Fourier  expansion  (7).  A  new  method  was 
proposed  by  Knab  et  al.  [17]  which  chooses  the  free  phase  factor  of  the 
CO-s  to  maximize  the  WF  norm  contribution  from  the  0-th  cell.  This 
method  gives  better  localization  than  the  original  one  implemented  in 
the  QP  program.  In  the  new  version  this  localization  is  built  in. 

•  The  loop  order  was  changed  (see  Fig  1).  In  the  new  version  the  outermost 
loop  is  for  the  WF  groups.  (  Non-redundant  WF  groups  are  defined 
on  the  same  way  like  the  non-redundant  AO  groups.  A  quadruplet  of 
cell  indices  connected  to  a  two-electron  integral  is  non-redundant  if  it  is 
not  possible  to  find  a  translation  or  a  trivial  permutation  of  cell  indices 
that  gives  an  other  member  in  the  non-redundant  set.)  The  integrals 
belonging  to  a  certain  WF  group  are  kept  in  the  memory  during  the 
transformation  and  written  out' when  the  trasformation  is  comleted  for 
the  particular  WF  group.  The  number  of  I/O  requests  was  reduced 
considerably  and  the  temporary  storage  is  also  decreased  removing  those 
work  files  where  the  partially  transformed  WF  integrals  were  stored.  In 
the  new  version  the  I/O  is  comparable  to  the  requirement  of  the  SCF 
calculation. 

•  Originally  the  different  steps  were  made  in  different  subroutines.  This 
way  especially  in  the  first  step  a  large  number  of  subroutine  calls  were 
carried  out  during  the  execution.  This  require  additional  computer  time. 
To  make  the  program  shorter  these  routines  were  built  in  to  the  code 
using  the  ’inline’  facility  of  the  Cray  Fortran  compiler. 
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•  To  calculate  the  self  energy  expression  we  do  not  need  the  total  two- 
electon  integral  list.  There  are  several  restrictions  for  the  bcind  indices 
in  eqn.  (3). 

—  We  need  only  those  integrals  where  the  first  index  is  equal  to  the 
index  of  the  band  to  be  corrected. 

—  In  the  first(second)  term  of  eqn.  (3)  the  summation  is  restricted  to 
two  occupied  and  one  virtual  band  index  and  to  two  virtual  and  one 
occupied  band  index,  respectively.  Further  only  the  different  pairs 
are  taken  from  the  occupied/virtual  subspace. 

When  (fo|uo)  and  {iv\ov)  type  integrals  are  transformed,  the  index  re¬ 
strictions  are  used  during  the  four  step  of  integral  calculation,  decreasing 
the  number  of  required  operations.  Here  {io\vo)  and  {iv\ov)  stand  for 
that  set  of  two-electron  integrals  where  the  indices  are  restricted  to  the 
occupied  (o)  and  virtual  (v)  bands  and  to  the  band  to  be  corrected  (i). 

•  Only  integrals  above  a  certain  threshold  value  are  stored  after  the  trans¬ 
formation  decreasing  the  required  peripheric  storage  capacity. 

•  A  restart  option  was  built  in  into  the  transformation  program.  It  is 
possible  to  restart  the  program  after  every  completed  WF  group. 


12 


Loop  for  HF  groups 
Loop  for  AO  groups 
Loop  for  / 

Loop  for  a 
Loop  for  (3 
Loop  for  7 
Loop  for  5 
1st  step 
Endloop  5 
2nd  step 
Endloop  7 


3rd  step 


Endloop  (3 
4th 

Endloop  a 
Endloop  I 
Endloop  WF 

Write  completed  WF  groups 
Endloop  AO 


Loop  for  oi  G  occ 

MI701) 

Endloop  Oi 
Loop  for  V\  G  virt 

(a/?|7Ui) 

Endloop  v\ 


Loop  for  vi  G  virt 
Loop  for  oj  G  occ 
(a/3|oiui) 
[a^\viOi) 

Endloop  Oi,Vi 


Loop  for  vi  G  virt 
Loop  for  oi  G  occ 
Loop  for  V2  G  virt 
(avi|oiV2) 
Endloop  V2tOi,vi 
Loop  for  oi  G  occ 
Loop  for  Vi  G  virt 
Loop  for  02  G  occ 
(aui|oi,U2) 
Endloop  02,vi,0i 


Loop  for  i  G  corrected  hands 
Loop  for  ui  G  virt 
Loop  for  oi  G  occ 
Loop  for  V2  G  virt  V2  > 
{ivi\oiV2) 

Endloop  V2,Oi,vi 
Loop  for  oj  G  occ 
Loop  for  ui  G  virt 
Loop  for  02  G  occ  02  >  Oi 
(:ui|oi,t;2) 

Endloop  i,02,ui,0i 


Figure  1:  Loop  struture  of  integral  transformation  program 


3.2  Applications 

3.2.1  PPP  calculations 

Poly(p-phenylene)  (PPP)  and  poly(p-phenylene  vinylene)  (PPV)  are  among 
the  most  promising  materials  for  optical  applications.  The  discovery  of  blue 
light  emission  from  PPP  and  fabrication  of  light  emission  diode  (LED)  made 
this  material  most  interesting  [19].  PPP  in  the  proper  structure  should  have 
a  band  gap  somewhat  larger  than  the  blue  frequency  regime  (due  to  the  most 
probable  occurrence  of  exciton  bands)  and  relatively  large  band  widths  of  the 
valence  and  conduction  bands  which  enable  sufficiently  high  mobility  of  the 
charge  carriers.  Several  theoretical  works  were  pulished  on  PPP.  Bredas  stud¬ 
ied  the  structure  of  PPP  oligomers  using  HF  and  semiempirical  methods  [20]. 
Ambrosch-Draxl  et  al.  [21]  carried  out  an  extended  local  density  functional 
investigation  of  ID  and  3D  PPP  structures  and  studied  the  possibilities  of  the 
tuning  of  gap  and  band  width  values. 

In  the  latest  Report  we  presented  the  first  neighbour  DZ  (dementi’s  DZ 
[24])  result  for  PPP.  Since  that  time  we  corrected  an  error  in  the  QP  program 
and  the  modified  results  are  presented  here.  We  used  the  same  geometry  as 
was  used  in  the  earlier  calculations  (Fig.  2). 


a=  1.404 
b=1.404 
c=l.S02 
d=  1.080 


Figure  2:  PPP  geometry.  Bond  distances  are  given  in  A,  bond  angles  are 
assumed  to  be  120°.  The  unit  cell  is  denoted  by  a  dashed  bracket 


In  the  first  case  we  used  three  interacting  cells  (first  neighbours’  interac¬ 
tions)  both  at  the  HF  and  QP  level  (Table  1)  .  The  valence  and  the  conduction 
bands  minima  cind  maxima  and  the  bandwidths  ((1)  are  presented.  The  last 
column  contains  the  absolute  value  of  the  QP  correction  and  the  arrows  (ff-IJ-) 
indicate  the  direction  of  the  shifts  of  the  bands.  The  last  row  contains  the 
HF  and  QP  gap  values  and  the  QP  correction.  The  bands  are  wide  ~4-6 
eV.  The  QP  gap  correction  is  as  large  as  3.36  eV.  The  QP  gap  vcJue  is  4.88 
eV  in  1st  neigbours’  and  5.63  eV  in  2nd  neighbours’  interactions  approxima¬ 
tion  Cctse  (optical  absorption  measurements  give  values  of  2.8  eV  [22]  cind  3.4 
eV  [23]  which  are,  due  to  the  presence  of  exciton  bands,  lower  bounds  of  the 
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fundamental  gap). 


HF 

QP 

Shift 

Conduction 

band 

max 

min 

5.09 

-0.61 

2.39 

-2.83 

$5.22 

$  2.70 
$  2.22 

Valence 

band 

max 

min 

-8.85 

-13.76 

$4.91 

-7.71 

-11.83 

$4.12 

$1.14 

$1.93 

GAP 

8.24 

4.88 

$  3.36 

Table  1:  PPP  DZ  calculation.  1st  neighbours’  interactions  approximation 
was  used  in  the  HF  and  QP  level.  The  values  are  in  eV 


In  our  second  calculation  we  used  five  interacting  cells  (second  neighbours’ 
interactions).  The  QP  gap  value  is  even  larger  than  in  the  previous  case.  The 
difference  is  0.75  eV,  this  comes  from  the  larger  HF  gap  (+0.46  eV)  and  the 
smaller  QP  correction  (-0.29  eV). 


HF  QP  Shift 

Conduction 

band 

max 

min 

5.57 

0.77 

$6.34 

4.93 

-1.34 

$6.27 

$  0.64 
$  2.11 

Valence 

band 

max 

min 

-7.93 

-12.45 

$4.52 

-6.97 

-10.78 

$3.81 

$  0.96 
$1.67 

GAP 

8.70 

5.63 

$  3.07 

Table  2:  PPP  DZ  calculation.  2nd  neighbours’  interactions  approximation 
was  used  in  the  HF  and  QP  level.  The  values  are  in  eV 

The  QP  corrections  at  the  lower  edge  of  the  conduction  band  and  at  the 
upper  edge  of  the  valence  band  differ  only  by  less  than  0.2  eV  in  the  first-  and 
second  neighbours’  interactions  cases,  respectively. 
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3.2.2  PPN  calculations 

Ladderlike  conjugated  hydrocarbons  play  very  important  role  among  conduct¬ 
ing  polymers.  A  Icirge  majority  of  these  polymers  are  insulators  in  spite  of 
their  mobile  x  electron  systems.  Poly(peri-naphtalene)  was  predicted  as  a 
good  intrinsic  conductor  fifteen  years  ago  [25].  The  valence  elfective  Hamil¬ 
tonian  (VEH)  calculations  predict  0.42-0.44  eV  as  band  gap  [26,  28],  minimal 
basis  HF  CO  calculation  gives  a  gap  vcdue  of  4.51  eV[27].  There  are  no  exper¬ 
imental  results  for  PPN  because  of  the  difficulties  in  the  synthesis.  Mullen  et 
al.  reported  an  extrapolated  0.92-0.98  eV  gap  value  using  optical  absorption 
measurment  results  of  PPN  oligomers  ((C'ioH4)„,n  =  2  —  5)  [29].  The  minimal 
basis  HF  result  overestimates  the  band  gap,  to  obtain  a  better  approximation 
we  need  a  better  basis  set  and  correlation  corrections. 

We  already  reported  the  minimal  basis  QP  calculation  for  the  PPN.  Now  we 
present  an  extended  basis  (dementi’s  DZ)  result  (Table  3).  First  neighbours’ 
interactions  were  taken  into  account.  We  used  a  PM3  optimized  geometry 
pulished  by  R.  Viruela-Martin  et  al.  [28]  and  showed  on  (Fig.  3). 

a=1.401 
b=1.382 
c=1.425 
d=1.415 
e=1.469 
h  =  1.10 

ab=121.00° 
be=119.3$° 
cd=119.64° 
ce=119.64  ° 
hb=119.44“ 


Figure  3:  PPN  geometry.  Bond  distances  are  given  in  A,  the  bond  angles  in 
degree.  The  unit  cell  is  denoted  by  a  dashed  bracket 

This  geometry  seems  to  be  reasonable  because  it  is  close  to  the  middle 
part  of  the  experimental  X-ray  diffraction  geometry  of  quaterrylene  [30]  that 
contains  four  naphtalene  units.  The  HF  gap  value  is  again  too  large,  but  with 
the  QP  correction  of  2.21  eV  the  gap  decreases  to  1.47  eV.  This  is  close  to  the 
predicted  experimental  gap,  the  difference  is  ~0.5  eV. 
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HF 

QP 

Shift 

Conduction 

band 

max 

min 

0.41 

-6.83 

$7.24 

-2.21 

-8.17 

$5.96 

$  2.62 
$  1.34 

Valence 

band 

max 

min 

-10.51 

-17.37 

$6.86 

-9.64 

-15.11 

$6.53 

$  0.87 
$  2.26 

GAP 

3.68 

1.47 

$  2.21 

Table  3:  PPN  DZ  calculation.  1st  neighbours’  interactions  approximation 
was  used  both  at  the  HF  and  QP  level.  The  values  are  in  eV 


3.2.3  Nucleotide  base  stacks  calculations 


DNA  plays  an  important  role  in  the  living  body  carrying  the  gemnetic  infor¬ 
mation.  Damages  in  DNA  can  result  in  cancer.  To  understand  the  processes 
that  cause  cancer  it  is  very  important  to  investigate  the  physical  structure  of 
this  molecule.  According  to  the  Wilkins- Watson-Crick  model,  DNA  forms  a 
double  helix  in  which  the  sugar-phosphate  backbones  are  outside  and  the  nu¬ 
cleotide  bases  (adenine  (A),  guanine  (G),  cytosine  (C),  thymine(T))  are  inside. 
The  sugar-phosphate  backbone  is  always  periodic.  Detailed  calculations  of  the 
backbone  [31]  have  shown,  that  the  valance  band  width  is  about  0.1  eV  and  the 
conduction  band  width  of  0.3  eV.  The  nucleotide  base  stacks  in  DNA  are  ape¬ 
riodic,  however  to  obtain  a  first  orientation  about  the  energy  level  distribution 
in  the  ba^e  stacks  in  DNA,  one  can  perform  calculations  of  the  energy  bands 
of  periodic  A  stack,  T  stack,  G  stack  and  C  stack.  These  so-called  periodic 
nucleotide  base  model  stacks  have  to  be  understood  such,  that  we  have  always 
the  same  nucleotide  base  repeated  in  the  same  relative  position  as  in  B  DNA 
(3.36  A  stacking  distance,  36°  rotation  [32]).  One  should  point  out,  however, 
that  on  one  hand  the  statistical  analysis  of  the  base  sequences  of  many  DNA 
fragments  has  shown  that  there  is  a  preference  in  native  DNA  to  have  the 
same  base  repeated  severah times  (in  a  few  cases  many  times  up  to  30).  On 
the  other  hand  all  the  so  called  homopolynucleotides  (a  sugar-phosphate  chain 
to  which  always  the  same  nucleotide  base  is  bound)  have  been  synthetized  in 
the  Laboratory  and  they  are  commercially  available. 

HF  level  calculation  have  already  been  reported  [33]  on  these  systems  using 
dementi’s  double  zeta  basis  set  [24]  .  Suhai  has  calculated  the  quasi-particle 
corrections  of  the  valence  and  conduction  bands  of  C  stack  [34].  Tables  4-7 
show  the  calculated  values  for  the  valence  and  conduction  bands  at  the  HF  and 
QP  level.  The  calculations  were  carried  out  using  the  dementi’s  double  zeta 
basis  set  [24]  and  in  the  first  neighbours’  interactions  approximation.  Again 
the  arrow  indicates  the  width  of  the  corresponding  band.  The  last  column 
contains  the  absolut  values  of  the  QP  corrections  at  the  band  edges,  Jj-  and  fT 
shows  their  direction. 

As  we  can  see  the  QP  corrections  are  very  important,  the  shift  of  the  va¬ 
lence  and  conduction  bands  is  1-1.5  eV,  and  the  decrease  in  the  gap  value 
is  2.5-3  eV.  The  smallest  QP  gap  value  is  8.19  eV,  it  shows  that  these  poly¬ 
mers  are  insulators.  Looking  at  the  band  widths  we  can  see  that  G  stack  has 
comparatively  the  largest  band  widths.  If  we  can  generate  free  charge  carri¬ 
ers  in  this  system  by  interaction  with  organic  electron  acceptors  or  donors  or 
electro  chemically,  the  G  stack  would  become  a  good  conductor. 

We  plan  to  extend  our  investigations,  using  higher  neighbours’  interactions 
approximation  and  larger  basis  sets  (polarization  functions).  Unfortunately 
the  inclusion  of  d-functions  means  5  additional  basis  function  for  each  non-H 
atom.  But  according  to  Qfzek  [35]  it  is  possible  to  select  a  proper  subset  of  d- 
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functions  which  provides  the  largest  part  of  d-functions  correlation  correction 
in  these  systems. 


HF 

QP 

Shift 

Conduction 

band 

max 

min 

3.21 

2.74 

ttO.47 

1.85 

1.51 

0.34 

$1.36 

$1.23 

Valence 

band 

max 

min 

-8.80 

-8.92 

0.12 

-7.13 

-7.26 

$0.13 

$1.67 
$  1.66 

GAP 

11.54 

8.64 

$  2.90 

Table  4:  C  stack,  DZ  calculation.  1st  neighbours’  interactions  approximation  B 

was  used  at  the  HF  and  QP  level.  The  values  are  in  eV 


HF 

QP 

Shift 

Conduction 

band 

max 

min 

3.73 

3.67 

$  0.06 

2.53 

2.20 

$  0.33 

$  1.20 
$  1.47 

VaJence 

band 

max 

min 

-8.38 

-8.60 

1 

$  0.22 

-7.32 

-7.38 

$  0.06 

$  1.06 
fr  2.22 

GAP 

-• 

12.05 

9.52 

$  2.53 

Table  5:  A  stack,  DZ  calculation.  1st  neighbours’  interactions  approximation 
was  used  at  the  HF  and  QP  level.  The  values  are  in  eV 


HF 

QP 

Shift 

Conduction 

max 

2.74 

1.38 

$  1.36 

band 

min 

2.53 

1.14 

$  0.24 

$  1.39 

Valence 

max 

-9.27 

-7.81 

$  1.46 

band 

min 

-9.98 

-8.37 

$  0.56 

$  1.03 

GAP 

11.80 

8.95 

$2.85 

Table  6:  T  stack,  DZ  calculation.  1st  neighbours’  interactions  approximation 
was  used  at  the  HF  and  QP  level.  The  values  are  in  eV 


HF 

QP 

Shift 

Conduction 

band 

max 

min 

4.11 

2.72 

3.26 

1.62 

$  1.64 

m 

Valence 

band 

max 

min 

-7.64 

-8.21 

$0.57 

-6.57 

-7.15 

$  0.58 

GAP 

10.36 

8.19 

$  2.17 

Table  7:  G  stack,  DZ  calculation.  1st  neighbours’  interactions  approximation 
was  used  at  the  HF  and  QP  level.  The  values  are  in  eV 
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3.2.4  Computational  times 


After  the  implementation  of  the  iterative  solution  of  the  inverse  Dyson  equation 
in  the  QP  program,  the  most  time  consuming  part  was  the  integral  transfor¬ 
mation.  For  PPP  in  DZ  basis  in  1st  neigbours’  interactions  approximation 
the  required  time  was  more  than  5000  second.  The  time  requirement  of  the 
different  parts  of  correlated  band  structure  calculation  is  shown  in  Table  8. 
PPPi  and  PPP2  means  the  1st  and  2nd  neighbours  PPP  calculation. 


Integral 

SCF 

Wannier 

Transf 

QP 

Total 

pppl 

2261 

1752 

65 

1699 

96 

5808 

PPP2 

11399 

2822 

130 

23455 

202 

38008 

PPN 

9565 

7073 

305 

8727 

535 

26205 

C  stack 

4302 

4450 

43 

4731 

274 

13800 

A  stack 

9224 

7123 

92 

11332 

574 

28345 

T  stack 

6425 

8301 

238 

7971 

423 

23358 

G  stack 

11477 

9548 

128 

16129 

727 

38009 

Table  8:  Computational  time  on  the  Cray  Y-MP8/8-128  (in  sec.) 

The  speed  up  factor  in  the  1st  neigbour  PPP  case  is  ~3.  We  can  see  that  in 
the  first  neighbours  case  the  integral  calculation  and  the  transformation  time 
is  comparable.  At  PPP  in  the  2nd  neighbours  case  the  transformation  time  is 
approximately  by  a  factor  of  two  larger,  than  the  integral  part. 
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